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Expectation-Maximization-Based Channel Estimation and 
Signal Detection for Wireless Communications Systems 

Field of the Invention 

The invention relates generally to communications systems, and more particularly 
to channel estimation and signal detection in wireless communications systems. 

Background of the Invention 

Orthogonal frequency division multiplexing (OFDM) is a signal modulation 
technique in which a transmitter divides a signal, and then transmits the divided 
signal over several subcarriers. The subcarriers are located on a frequency axis at 
regular intervals. With the OFDM technique, in contrast with conventional serial 
communication techniques, the transmitted signal is divided into N streams, and 
the N stream are then transmitted in parallel over N subcarriers, each having a 
separate carrier frequency. The OFDM technique transmits the signal reliably and 
efficiently at a high data rate. 

The subcarriers are made "orthogonal" by appropriately selecting the spacing of 
the frequencies in the frequency band. Therefore, spectral overlapping among 
subcarriers is allowed because the orthogonality ensures that the receiver can 
separate the OFDM subcarriers. With OFDM, a better spectral efficiency is 
achieved than by using a simple frequency division multiplexing technique. OFDM 
is more robust to data loss due to multipath fading when compared with a single 
carrier because OFDM has an increased symbol period for the same aggregate data 
rate. 
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In addition, inter-symbol interference (ISI) in OFDM transmissions can be 
prevented by inserting a guard interval before each transmitted block of symbols. 
Moreover, OFDM is robust to frequency selective fading because each sub-channel 

5 occupies a relatively narrow frequency band, where the characteristic of the 

channel frequency is relatively flat. Thus, OFDM is used by many communication 
systems, including digital audio and video broadcasting (DAB, DVB), and high- 
speed digital subscriber line (DSL) modems over a twisted pair of wires. OFDM 
can also be used in wireless local area networks (WLANs), and fixed broadband 

W wireless communication networks. 

£ However, it is not possible to make reliable data decisions in OFDM systems 
tO unless a good channel estimate is available. Thus, an efficient and accurate channel 
n estimation method is necessary to coherently demodulate received data. Although 
Q differential detection could be used to detect the transmitted signals in the absence 
S of channel information, this would result in about a 3dB loss in SNR compared to 
m coherent detection. 

A number of channel estimation techniques are known in the prior art. In most of 
20 those, the channel estimates are continuously updated by transmitting pilot 

symbols using specified time-frequency lattices. One class of such pilot assisted 
estimation processes adopt an interpolation technique with fixed one or two- 
dimensional parameters to estimate the frequency domain channel impulse 
response (CIR). Channel estimates are obtained at lattices assigned to the pilot 
25 tones, see Jae Kyoung Moon and Song In Choi, "Performance of channel 
estimation methods for OFDM systems in multipath fading channels," IEEE 
Transactions on Consumer Electronics, Vol. 46, No. 1, Feb. 2000, pp. 161-170, P. 
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Hoeher, S. Kaiser, and P. Robertson, "Two-dimensional pilot-symbol-sided 
channel estimation by Wiener filtering," Proceedings of 1997 IEEE International 
Conference on Acoustics, Speech and Signal Processing, ICASSP-97, vol. 3, pp. 
1845-1848, and F. Said and A. H. Aghvami, "Linear two dimensional pilot assisted 
5 channel estimation for OFDM systems," 6 th IEEE Conference on 

Telecommunications, 1998, pp. 32-36. Linear, spline and Gaussian filters can all be 
used with these methods. 



Another class of methods adopt known channel statistics and channel estimates at 
ffi pilot symbols to estimate the CIR in the sense of minimum mean square error 
9 (MMSE), see Ye Li, Leonard J. Cimini, Jr., and Nelson R. Sollenberger, "Robust 
%l channel estimation for OFDM systems with rapid dispersive fading channels," 
IEEE Transactions on Communication, Vol. 46, No. 7, My 1998, pp. 902-915, J.- 
J. van de Beek, O. Edfors, M. Sandell, S. K. Wilson, and P. O. Borjesson, "On 
\S channel estimation in OFDM systems," IEEE Vehicular Technology Conference, 

m 1995, Vol. 2, pp. 815 -819, and O. Edfors, M. Sandell, S. K. Wilson , J.-J. van de 

□I 

fjj Beek, and P. O. Borjesson, "OFDM channel estimation by singular value 

decomposition," IEEE Transactions on Communications, Vol. 46 No. 7, July 1998, 
pp. 931 -939. Shortcomings of those processes include a large error floor that may 

20 be incurred by a mismatch between the estimated and real CIRs, and difficulty in 
obtaining the correlation function of the channel impulse response. 



Therefore, there is a need for an efficient method for estimating channels so that 
received symbols can be coherently detected and demodulated. 

25 
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Summary of the Invention 

The invention provides an expectation-maximization (EM) based system and 
method for efficiently estimating channel states and for coherently detecting 
5 received symbols via channels of a wireless communications system. A maximum 
likelihood estimate of the channel impulse response (CIR) is obtained by using the 
EM-based method to estimate the channel mean and covariance, i.e., the channel 
statistics. The invention improves signal detection, and the accuracy of the channel 
estimate by making use of the iterative EM process. 

sash 

5 

p More particularly the system and method detects symbols of a modulated signal 

VI 

j*& received via a plurality of channels of a wireless communications system. The 

y§ method updates the channel estimate and detects the received symbol iteratively 

q until the estimates converge. The estimate of the symbol is first determined by 

S M 

fl maximizing the expectation of the log likelihood function, and then quantized 
JJ: according to the signal constellation. 

s 

The quantized estimate of the symbol is compared with the previous estimate of 
the symbol to determine if the estimates have converged. Otherwise, the quantized 
20 estimate is made as the input for the next iteration, and the estimating, maximizing, 
quantizing, and comparing are repeated until the estimates converge. 

Brief Description of the Drawings 

25 Figure 1 is a block diagram of a wireless communications system according to the 
invention; 
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Figure 2 is a flow diagram of a general method for detecting transmitted symbols 
in the system of Figure 1 ; 

Figure 3 is a flow diagram of a simplified method for detecting symbols in the 
5 system of Figure 1; 

Figure 4 is a flow diagram of a basic method for detecting symbols in a MPSK 
modulated system; 

W Figure 5 is a graph comparing bit error rates for different channel estimation 
in methods as a function of SNR; 

di Figure 6 is a graph comparing mean square error for different channel estimation 
n methods as a function of SNR; and 

ft ' 

SI Figure 7 is a graph comparing iterations for signal detection as a function of SNR. 

fPJ 

Detailed Description of the Preferred Embodiment 

20 Expectation-Maximization 

Our invention uses an expectation-maximization (EM) method to estimate 
statistical characteristics (parameters) of a wireless communications channel, and 
to detect signals in the estimated channel. 

25 

The EM method according to the invention includes two iterative steps; an 
expectation step and a maximization step. The expectation step is performed with 
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respect to unknown underlying channel parameters, using a previous estimate of 
the parameters, conditioned upon incomplete observations. The maximization step 
then provides an improved estimate that maximizes the expectation of a log 
likelihood function defined over complete data, conditioned on the previous 

5 estimate and the next estimate. These two steps are iterated until the estimated 
values converge. See A. P. Dempster, N. M. Laird and D. B. Rubin, "Maximun 
likelihood estimation from incomplete data," Journal of the Royal Statistical 
Society (B), Vol. 39, No. 1, 1977, pp. 1-38, and T. K. Moon, "The expectation- 
maximization algorithm," IEEE Signal Processing Magazine, Vol. 13 No. 6, Nov. 

tO 1996, pp. 47 -60 for a general description of EM-based methods. 

We use EM to estimate channel states and to recover transmitted symbols of an 
% £ OFDM system over a frequency selective channel with additive white Gaussian 

noise (AWGN). After each iteration, we obtain a more accurate channel estimate 
%% for detecting transmitted symbols. 

IS ii 

U OFDM System Structure 
Transmitter 

20 

Figure 1 shows a baseband equivalent representation of an OFDM system 100 that 
uses the invention. In a transmitter of the OFDM system 100, input binary data 
symbols X 101 are fed into a serial to parallel (S/P) converter 1 10 to produce 
parallel data streams. Each parallel data stream 111 is then modulated (Mod) 120 
25 to corresponding sub-carriers 121 using, for example, MPSK or MQAM 

modulation techniques. Modulation schemes can vary on the sub-carriers in order 
to achieve a maximum capacity or a minimum bit error rate (BER) under some 
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constraints. We describe the invention in terms of fixed modulation, e.g., QPSK or 
16QAM, on all of the sub-carriers 121. However, it should be understood that the 
invention could be applied to any modulation scheme. 

5 The modulated symbols, represented by complex variables X (0), • • • , X (M - 1) , are 
transformed by an inverse fast Fourier transform (IFFT) 130. The transformed 
symbols are denoted as x(0),- • *(M-1) . In order to avoid inter-symbol interference 
(ISI), a cyclic prefix (CP) is added 140 in front of each symbol. The CP replicates 
M the end part of the previous IFFT output symbol. The length of the CP is longer 
III than the channel length L. The parallel data are converted (P/S) 150 to a serial data 
M stream 105, which is transmitted over a frequency selective fading channel 102. 

f " Receiver 

jj In a receiver of the system 1 00, received symbols y (0), — , y(M - 1) 1 03 corrupted by 
m multipara fading and AWGN, are converted (S/P) 160 to a form F (0), • • • , Y(M - 1) , 
after the cyclic prefix is discarded 165. A fast Fourier transform (FFT) 170 is 
applied and the signal is demodulated 180 and converted 190 to serial form as 
output symbols X 104 that correspond to the input symbols 101. 

20 

Channel Model 

In our channel model, the variables H,h,X_,Y,N denote the vectors of frequency- 
domain OR, time-domain OR, transmitted symbol, received symbol and white 
25 Gaussian noise, respectively. Our multipath time-invariant fading channel 102 is 
described by 
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1-1 

yik)=Y/hx(k-l)+Ti(k), 0<k<M-l f (1) 

1=0 

where the time-domain CIR, V s (0</<L-l) are independent complex-valued 
Rayleigh fading random variables, and n(k) 's (0 < k < M - 1) are independent 
complex-valued Gaussian random variables with a zero mean and a variance o 2 
5 for both real and imaginary components. The variable L is the length of the time- 
domain CIR. 



With the added CP in each OFDM symbol, there is no ISI. We consider only one 
OFDM symbol with M sub-carriers in analyzing the model of our system. After 
discarding the cyclic prefix and performing the FFT 170 at the receiver, we obtain 
the received symbol in the frequency domain as: 

i M-l _ J2n bn 

Y(m)=-j=^j(k)e \ (2) 



b=0 



O Substituting (1) into (2), we have 

15 Y(m)=X(m)H(m)+N(m), 0<m<M-\, (3) 

where H(m) is the frequency response of the channel 102 at a sub-carrier m, which 
can be expressed as 



r_] ._ ml 

1 ~j2tz — 



tf(m)=£V M ' 0<m<M-l. (4) 



/=0 



20 A set of the transformed noise variables { N(m), 0 < m < M - 1 } , which can be 
written as 

h M-l _ 2ir— 

N(m)=^=J/i(k)e J M , 0<m<M-h (5) 
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are i.i.d. complex-valued Gaussian variables that have the same distribution as 
{ n(k) }, i.e., with zero mean and variance o 2 . Note that inter-carrier interference 
(ICI) is eliminated at the output of the FFT 170 because of the orthogonality 
between the sub-carriers. 

5 

EM-Based Channel Estimation and Signal Detection Method 

Our objective is to detect the transmitted symbol X(m), 0 < m < M - 1 105 from the 
received or observed symbol Y(m), 0 < m < M -1 103. In order to reduce the bit error 

W rate caused by uncertainty in the channel 102, we apply our EM-based method to 

p 

jjfl take an average over the unknown CIR, assuming that the probability density 

1=2 function (pdf) of the CIR response is known at the receiver. 

JU As previously stated, each transmitted symbol X(m) 101 is modulated, for 

jl example, by QPSK or 16QAM. We use the notations *, T and H to denote the 

Ji; complex conjugate, transpose and complex conjugate transpose (Hermite), 

m respectively, in our analysis. 



Thus, we have h = [/t, , • • ■ , h, f , X = [X (0), • • • , X (M - l)f , Y = [Y(0), • • • , Y(M - l)f , 
20 iV = [iV(0), --,iV(M-l)] 7 ', H_ = Wh , and W is a MxL matrix: 



W = 



- j2x — 

1 e M 



1 e M 



-jlTT 



(Af-1) 



M 



.^(M-JKL-l) 
9 M 



(6) 
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We use the notation X = diag(X) to denote a MxM matrix with values X at 
diagonal terms of the matrix, and zeros elsewhere. 

We assume there is no ISI between two successive OFDM symbols because the 
5 cyclic prefix is longer than the channel length L. Thus, we need only to consider 
one OFDM symbol for our analysis, and therefore omit the symbol index, and 
express the channel model by 

Y = XWh + N. (7) 

W Our EM-based method obtains an estimate of the transmitted symbol X that 

maximizes likelihood p(Y I X) by averaging the logarithm of a likelihood function 
Ci! PQL> h 1 20 over th e unknown channel parameters h . 



v. 



The "incomplete" and "complete" symbols are (7) and (Y, h) , respectively. At each 
jjf iterative step p(p =0, X % • ••) , the EM method for estimating the transmitted 
|j symbol X from the received symbol Y includes the following two steps: 
E-step: 

Q(X\K p ) = E{\og p(Y, h\x)\Y,X p }, and (8) 
M-step: 

20 K p+i =argmax0(X | X p ). (9) 

During the (p + l) s ' iteration, the E-step determines the expected log likelihood 
function Q(X \X p ), given Y and X p which are the estimates obtained during the 
previous p ,h iteration. 

25 
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During the same iteration, the M-step determines the transmitted symbol X p+l that 
maximizes Q(X | X p ) for the next step, given the previous estimate of the symbol 



X_ p . 



Equation (8) of the E-step can be rewritten as 

QQL | X p ) = Jlogp(L h | X)p(h | F, X p )dh , (10) 
where the log likelihood function can be expressed as 

log P(L h\X) = log p(Y \KX) + log p(h 1 2D. (1 1) 



The conditional pdf 1 7, X p ) used in equation (10) takes the conditional 
expectation over the unknown channel parameters h . We assume that the CIR h 
and X are independent of each other. This is a reasonable assumption because, in 
general, the CIR does not depend on the transmitted symbol. Thus, for the purpose 
of maximizing equation (9), the Q function of equation (10) can be replaced by 

Q(K | X p ) = JlogpO: | h, X)p(h I Y, X p )dh . (12) 

The logarithm of the likelihood function p{h \ Y, X_ p ) can be determined by 
. p(Y\h,XJp(h\X p ) p(Y\h,X p )p(h) 



where we use the assumption that h and X_ p are independent of each other. Thus, 
equation (12) can be reduced to 

Q(X | X p ) = jlogpd | h, X)p(Y | h, X p )p(h)dh . (14) 
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Because the likelihood p(Y \ X p ) does not depend on X , it can be discarded in the 
last expression. 

We now determine the exact expression of Q(X_ | X p ) for a fading channel with 
AWGN. The conditional pdf s p(Y \ h, X) and p(Y \h,X p ) take the form 

Pd\ A,X) = (W)- M expH|7-X\^|72c7 2 } j ^ 

Pd | b K p ) =(W)- M expf ^-X^f/V} 

where o 2 is the variance of both the real and imaginary components of the 
complex-valued Gaussian noise. The pdf p(h)is given by 

pQi) = (^Idetl)" 1 ex^ih-EfeYx-'ih-Efe})} 

where E{h} and S are the mean and covariance matrix of the CIR h . 



By omitting the normalization constants, equation (12) can be expressed as 
e(2C | K p ) = -J ||y - xw/*|| I y, x p )dh , (18) 

and p(h \Y,X_ p ) can be represented as 

p(h\LX p ) = Kcx 1? \-^h-h p ) H t- p \h-h p )l (19) 
where K is a normalized constant. The values h p and £ p are called the estimated 
posterior mean and posterior covariance matrix at the p ,h iteration, respectively. 



These can be expressed as 

h p =± p (W H X H p Y/(7 2 + I.- 1 Efe} y md (20) 

t p = (W h X h p X p W/ct 2 . (21) 
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Maximization of equation (18) is equivalent to minimizing a distance in an 
integral, i.e., 

arffDS30X\X p )^ (22) 

.^c x 

This minimization can be further simplified to 

mzxE^ H F+F H h-h H Gh \Y,X p \ (23) 

where 

F = W W X // F, and (24) 
G = W f, X fl XW. (25) 



M Because the distribution of the random vector Y, given h and X , is Gaussian 
with a mean ft p and covariance matrix t p , it is possible to determine 



Moreover, all entries of the matrix G are given in terms of the signal energies, i.e., 
2 ,---,||X(M -l)f. Thus, we determine the third part of equation (19) by 

M 



E&GhlLX^CfWXdf, (27) 

1=1 

where values of a channel update coefficients matrix C?, 0 < / < M - 1 are real 
values dependent on h p mdt p . These values are obtained according to the 
20 following equation: 

Cf = (28) 



m=Q n=0 
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In order to determine Q(X \ X p ) completely, we rewrite equation (24) as follows 
F = W H YX\ (29) 

Thus, maximizing Q(X_ \X_ p ) is the same as 

max |a" W h YX_ * + X T YW h p - £ C, 2 |X(f)| 2 J (30) 

Next, we maximize equation (30), in terms of the transmitted symbol X , by 
solving 

£ p+1 =argmaxQ(X |X p ) = C- , ^W w Y) r , (31) 

where C = diag(C 0 C M _ X ) . After quantizing X_ p+ , , we obtain the next (p + 1)* 
estimate 

X_ P +\ = Quantization^ p+l )• (32) 

At each iteration, the updated estimate of CIR h p is obtained automatically as a 

by-product. The next estimate of the symbol is then compared with the previous 
estimate of the symbol to determine if the estimates have converged. On 
convergence, estimation of the transmitted symbol is complete, otherwise reiterate. 

General Method EM-Based Channel Estimation and Signal Detection 

Figure 2 shows the general EM-based iterative signal detection method according 
to the invention. 
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For step p = 0, an initial estimate of the transmitted symbol, X_ 0 , is obtained 200 
for the received signal Y 201 by using the received signal Y 201 and the channel 
estimate 202 from a pilot symbol inserted in an OFDM frame, or from the 
previous estimated channel information under the assumption that the channel 
varies slowly. 

Then for each iteration p, first determine 210 a posterior covariance matrix L p of 
the channel using the EFT matrix W, the previous p th symbol estimate X p , which is 
X_ 0 initially, the channel convergence matrix 2T 1 , and the Gaussian noise variance 
o 2 as: 

± p =(W H XjX ; ,W/<7 2 +E- , r 1 . 
In step 220, determine the posterior mean h p of the channel impulse response 
using the received data, i.e. the observed symbol Y. 

In step 230, determine the channel update coefficients matrix C for recovering 
estimated symbols by evaluating 

M-1M-1 f 2, T (m " n)/ / \ 

C ( 2 =XI> M $ P (m,n) + h;(m)h p (n)). 

m=0 n=0 

Then, in step 240, determine the estimate for the (p + l) st iteration X p+1 by 
applying the coefficient matrix C to the posterior mean h p of the CIR, the FFT 
matrix W, and the received symbol Y according to 
l P+ i =C- 1 (h H p W H Y) T . 
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Step 250 quantizes the next estimate of the symbol according to the signal 
constellation 



In step 260, compare the next (p + l) st estimate X_ p+1 with the previous p st estimate 



the estimate of the transmitted symbol has converged 270, and the iterations stop. 
Otherwise reiterate beginning at step 210 until X_ p+1 ~ X_ p , that is, until result of the 

next iteration is substantially equal to the result of the previous iteration. 
Simplified EM Method 

Thus far, we have assumed that the length of the channel L, is known. In a real 
situation, however, L is generally not known. In such a case, we need to perform 
channel-order detection together with parameter estimation. 

Alternatively, we can use some upper bound for L, which may be easier to obtain 
than trying to estimate an exact value of L. For the OFDM system 100, we can set 
L equal to the length of the cyclic prefix because, as we stated before, the cyclic 
prefix is at least longer than the channel delay spread in order to eliminate ISI. 



Kp+i = Quantization 




X„ . If an absolute difference between the estimates is smaller than a 



p 



predetermined threshold, 
\X p+l -X p \<T, 
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In addition, the memE{h} and covariance matrix £ of the time-domain CIR also 
need to be known. In the practical situation, these channel statistics are difficult to 
obtain. 

As is known from the general convergence property of the EM method, there is no 
guarantee that the iterative method converges to a global maximum. For a 
likelihood function with multiple local maxima, the convergence point may be any 
one of these local maxima, depending on the initial estimated symbol X 0 . 

Therefore, we use pilot symbols distributed at certain locations in the OFDM time- 
frequency lattice to obtain appropriate initial values of the CIR and X 0 which are 
most likely to lead us to the true maximum. 

So far, we have assumed that the channel statistics, i.e., the mean and covariance 
matrix, are known at the receiver. However, as we stated above, the channel 
statistics are difficult to obtain in practical situations. Fortunately, when we 
examine equations (26) and (27), we find that when the variance a 2 is small, i.e., 
the SNR is high, the contribution of L is so small that we can eliminate it and still 
obtain a similar performance. 

Therefore, as shown in Figure 3, we can simplify the general EM-based iterative 
signal detection method as follows. 

As before in the initial step p = 0, an initial estimate of the transmitted symbol, X 0 , 
is determined 300. The initial estimate can be obtained by using the received signal 
Y and the channel estimate from a pilot symbol inserted in the OFDM frame, or 
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from the some previous estimated channel information under the assumption that 
the channel varies slowly. 

Then for each iteration p, first determine 310 a posterior covariance matrix L p of 

5 the channel using the FFT matrix W, the previous p th symbol estimate X p , which is 
X 0 initially, and the Gaussian noise variance o 2 as: 

Note that in this simplified method the channel covariance matrix I" 1 is not used. 

f§ In step 320, determine the posterior mean h p of the channel impulse response 
H h p =± p (W H X H p Y/a 2 ), 

M using the received data, i.e., the observed symbol Y. Note that in the simplified 

« method the CIR E^h} is also not used. 

|| In step 330, determine the channel update coefficients matrix C for recovering 

O 

f|| estimated symbols by evaluating 

C-=Y%e " $ p (m,n) + h;(m)h p (n)). 

m-0 n-0 

Then, in step 340, determine the estimate of the transmitted symbol Z p+l for the (p 

20 + l) s iteration by applying the coefficient matrix C to the posterior mean h p of the 
CIR, the FFT matrix W, and the received signal Y according to 
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Step 350 quantizes the next estimate of the transmitted symbol according to the 
signal constellation 

2Lp+i = Quantization(x_ p+l ). 

5 In step 360, compare the next (p + l) st estimate X_ p+l with the previous ^'estimate 
X p . If an absolute difference between the estimates is smaller than a 

predetermined threshold, 
\X p+l -X p \<T, 

jj; the estimate of the transmitted symbol has converged 370, and the iterations stop. 
I§ Otherwise reiterate beginning at step 310 until X_ p+1 »X P , that is, result of the next 

J; iteration is substantially equal to the result of the previous iteration. 

m 

k Basic EM-Based Method for MPSK Modulated Signal 

j§ For a MPSK modulated signal, i.e., ||X(m)|| 2 = A for all m , where A is a positive 

m constant equivalent to the signal energy, the estimation can be performed by using 
only phase information. 

Thus, as shown in Figure 4, we can further simplify the above method as follows. 

20 

As before in the initial step p = 0, an initial estimate of the transmitted symbol, 
X 0 , is determined 410. The initial estimate can be obtained by using the received 
signal Y and the channel estimate from a pilot symbol inserted in the OFDM frame, 
or from the previous estimated channel information under the assumption that the 
25 channel varies slowly. 
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Then for each iteration p, determine the estimate for the (p + l) st iteration X p+1 by 

using only the FFT matrix W, the received signal Y, and the previous p th symbol 
estimate X p according to 

X p+l =(YX p WN H Y) T . 

Step 430 quantizes the next estimate of the symbol according to the signal 
constellation 

X_ p+l = Quantization^ p+} \ 

In step 440, compare the next (p + l) st estimate with the previous p st estimate 

X . If the absolute difference between the estimates is smaller than a 

p_ 

predetermined threshold, 
\X p+l -X p \<T, 

the estimate of the transmitted symbol has converged 450, and the iterations stop. 
Otherwise reiterate beginning at step 410 until X p+l ~ X p , that is, result of the next 

iteration is substantially equal to the result of the previous iteration. 

In this basic form, we need only to perform addition and multiplication operations, 
and WW H can be determined ahead of time. 
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Effect of the Invention 

We can construct an OFDM model system to demonstrate the validity and 
effectiveness of the EM-based signal estimation method according to the invention. 
5 The entire channel bandwidth of 400kHz is divided into 64 subcarriers. 

The symbol duration of each orthogonal subcarrier is 160 jus . An additional 20 jus 
cyclic prefix added to provide protection from ISI and ICI due to channel delay 
spread. Thus, the total OFDM frame length is 7>180 jus , and the subcarrier 
If) symbol rate is 5.56 kilobaud. 

In this model of the system, the modulation scheme is QPSK. One OFDM frame 



out of eight OFDM frames includes pilot symbols, and eight pilot symbols are 
inserted into such a frame. The system can transmit data at 700 kbits per second. 



IS The maximum Doppler frequencies are 55.6 Hz and 277.8 Hz, which make// 7; 
jlj 0.01 and 0.05, respectively. The channel impulse response is 



where C - Y^_e i s a normalization constant and a k 0 <k<l are independent 



complex-valued Rayleigh distributed random variables with unit energy, which 
20 vary in time according to the Doppler frequency. This is a conventional 
exponential decay multipath channel model. 

Figure 5 shows the BER performance of EM-based OFDM signal estimation 
method with the above two different Doppler frequencies, and Figure 6 shows the 



h{n)=-^e- kl2 a k 8{n-k) y 



(33) 
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corresponding MSE (Mean Square Error) of the estimate. In both graphs the 
performance is shown as a function of the SNR E b /N 0 . 

In prior art channel estimation methods, for those OFDM frames containing pilot 
5 symbols; the estimate of CIR is obtained by using these eight equally spaced pilot 
symbols. For those OFDM frames without pilot symbols, the conventional 
estimation of CIR comes from the channel estimate of the previous OFDM frame. 

In the EM-based method according to the invention, the channel estimate of the 
10 previous symbol is used for obtaining the initial value for the next symbol when 
p there are no pilot symbols in the frame. 

Q From Figures 5 and 6, it can be seen that our EM-based method achieves 

substantially the same performance as the ideal case in terms of BER, where the 
| J channel characteristics are completely known when/ d T s = 0.01 , i.e., the channel 
[U characteristics do not change very fast. Furthermore, the MSE of our EM-based 
O method converges to the Cramer-Rao Lower Bound (CRLB) when E b /N 0 becomes 

large. 

20 The performance gain from the initial estimate is considerably large, especially 
when/j T s = 0.05. In addition, the performance degradation is substantially small 
when we use the simplified EM-based method that does not use the channel 
statistics. Thus, our method is robust. 

25 In Figure 7, we show the number of iterations required for the estimates X_ p to 
converge versus E b /N 0 at the receiver input. We see that the number of necessary 
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iterations is relatively small in a broad range of E b /N 0 , whether the channel varies 

slowly, or not at all. Fast channel variation causes only a very small increase in the 
number of iterations required to converge. This demonstrates that our method can 
achieve a substantial performance improvement with only a modest increase in the 
computational complexity. 

We have described a novel EM-based channel estimation and signal detection 
method, its simplified form for OFDM systems, and a basic form for MPSK 
systems. The method is efficient because the number of iterations needed to 
converge is small, e.g., three or less, and the complexity of the operations required 
in each iteration is low. By making use of pilot symbols to obtain the initial 
estimate, the method can achieve a near-optimal estimate after a small number of 
iterations when the channel changes slowly. The MSE of the channel estimation 
approaches CRLB when E b /N 0 > XOdB . 

The method according to the invention can reduce the number of pilot symbols 
inserted in the time-frequency grid of an OFDM system over slow time variation 
channels. Thus, the spectral efficiency is improved. The corresponding simplified 
method does not need any knowledge of the channel statistics without degrading 
performance. 
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The method according to the invention can also improve channel estimation 
accuracy obtained from any other method by iterative steps using that channel 
estimate as the initial estimate. Our method can also be used to estimate the 
channels in multi-input/multi-output (MIMO) OFDM systems. 

Although the invention has been described by way of examples of preferred 
embodiments, it is to be understood that various other adaptations and 
modifications may be made within the spirit and scope of the invention. Therefore, 
it is the object of the appended claims to cover all such variations and 
modifications as come within the true spirit and scope of the invention. 
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